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ABSTRACT 

As is well known, a varying effective gravitational "constant" is one of the common features of most 
modified gravity theories. Of course, as a modified gravity theory, f(T) theory is not an exception. 
Noting that the observational constraint on the varying gravitational "constant" is very tight, in 
the present work we try to constrain f(T) theories with the varying gravitational "constant". We 
find that the allowed model parameter n or /3 has been significantly shrunk to a very narrow range 
around zero. The results are very impressive. 
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I. INTRODUCTION 

Motivated by the well-known large number hypothesis proposed in 1937 [l|, the varying fundamental 
"constants" remain as one of the unfading subjects for decades. In fact, the gravitational constant G was 
the first one whose constancy was questioned by Dirac [l[ . 

As a pure gravitational phenomenon, the variation of the gravitational "constant" does not affect the 
local physics (e.g. the atomic transitions or the nuclear physics), and hence most constraints on it are 
obtained from systems in which gravity is non- negligible, such as the motion of the bodies of solar system 
and astrophysical systems [2| • Following [2| , here we briefly review the current observational constraints 
on the varying gravitational "constant" . The first type of constraints arises from solar system. The latest 
analysis of the lunar laser ranging experiment Q gives the constraint 

\G/G\ < 1.3 x 10- 12 yr -1 , (1) 

where a dot denotes the derivative with respect to the time t. In [4], the combination of Mariner 10, 
Mercury and Venus ranging data gives 

\G/G\ < 2 x KT 12 yr -1 - (2) 

In [f| , the ranging data from Viking landers on Mars leads to the constraint 

\G/G\ < 6 x HT^yr- 1 . (3) 

The second type of constraints arises from pulsar timing. In |6[, the constraints from PSR B1913+16 
and PSR B1855+09 are given by 

\G/G\ <9x 10- 12 yr _1 , (4) 

and 

\G/G\ < 2.7 x 10 _11 yT _1 , (5) 
respectively. On the other hand, the constraint from PSR J0437-4715 reads Q 

\G/G\ <2.3x HT 11 ^" 1 - (6) 
The third type is the stellar constraints. In Q, the ages of globular clusters give the constraint 

\G/G\ <3.5x HT 11 ^" 1 - (7) 
The constraint from helioseismology is given by [9] 

\G/G\ <1.6x lO-^yr" 1 . (8) 
The seismology of white dwarf G117-B15A [l(| gives the constraint 

\G/G\ <4.1 x HT 11 ^" 1 - (9) 
In [ll[ , the constraint from the cooling of white dwarfs reads 

|G/G|<2xlO _11 yr~ 1 - (10) 
In [l2j , the light curves of supernovae give the constraint 

\G/G\ <4.8x lO -12 ^- 1 . (11) 

Finally, one can also obtain the cosmological constraints on the varying G from cosmic microwave back- 
ground (CMB) and big bang nucleosynthesis (BBN). However, this type of cosmological constraints on 
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the varying G heavily relies on not only the assumption of the whole history of G(t) but also many other 
complicated factors. Therefore, it is very difficult to state a definitive constraint [2j. So, here we do not 
consider this type of the cosmological constraints on the varying gravitational "constant" . 

Obviously, the currently tightest observational constraint on the varying gravitational "constant" is 
the one given in Eq. ([1]), namely, \G/G\ < 1.3 x 10 _12 yr _1 . So, in the present work we only use this 
tightest observational constraint given in Eq. ([1]) to constrain f(T) theories. 

In fact, it is well known that a varying effective gravitational "constant" is one of the common features 
of many modified gravity theories [2|, for examples, f(R) theory, scalar-tensor theory (including Brans- 
Dicke theory and Galileon theory), braneworld scenarios (such as DGP, RSI and RSII), f(Q) theory (Q is 
the Gauss-Bonett term), Horava-Lifshitz theory, MOND and TeVeS theories (see e.g. (13l - [l7| for reviews). 
Recently, a new modified gravity theory, namely the so-called f(T) theory, attracted much attention in 
the community, where T is the torsion scalar. Similar to many modified gravity theories, the effective 
gravitational "constant" is also varying in f(T) theory. Therefore, it is of interest to constrain f(T) 
theory with the varying gravitational "constant" . This is what we need to do in the present work. 

In Sec. HH we briefly review the key points of f{T) theory. In Sec. MI Al in order to compare with the 
observational constraint on the varying gravitational "constant" , we give the corresponding formula of 
|G e ff/G c ff| for a general f(T) theory. In Sees. IIII Bl and ITlI G[ we constrain two concrete f(T) theories 
with the varying gravitational "constant" , respectively. Finally, a brief conclusion is given in Sec. IIV1 



II. A BRIEF REVIEW OF f(T) THEORY 

f(T) theory is a generalization of the teleparallel gravity originally proposed by Einstein 0>E|- I n 
teleparallel gravity, the Weitzenbock connection is used, rather than the Levi-Civita connection which is 
used in general relativity. Following [5q[[21|, here we briefly review the key points of f(T) theory. The 
orthonormal tetrad components ej(x M ) relate to the metric through 

9»u = Vij^ei , (12) 

where Latin i, j are indices running over 0, 1, 2, 3 for the tangent space of the manifold, and Greek fi, v 
are the coordinate indices on the manifold, also running over 0, 1, 2, 3. In f(T) theory, the gravitational 
action is given by 

St = j^I d4x]e] [T+/(T)] ' (13) 

where « 2 = 8ttGn (Gn is the Newton constant), and |e| — det (e^) — \J~g ■ It is worth noting that in 
the literature, T + f(T) in Eq. (fT3")) could be instead replaced by f(T), and hence one should be aware 
of the correspondence between these two formalisms. The torsion scalar T is defined by 



T = S i»> T p (14) 



where 



TV = -< (Vt " , (15) 

K^ p = -~ (T»% - T v % - T/") , (16) 

Sr = \ (K»\ + <5£T% - S^T 9 ^) . (17) 

In this work, we consider a spatially flat Fricdmann-Robertson- Walker (FRW) universe. In this case, one 
can easily find that [2(J Hl| 

T = -QH 2 , (18) 

where H = a/a is the Hubble parameter; a — (1 + z)^ 1 is the scale factor (we have set ao=l, whereas 
the subscript "0" indicates the present value of the corresponding quantity); z is the redshift. Therefore, 
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one can use T and H interchangeably. The modified Friedmann equation and Raychaudhuri equation are 
given by [HHH 

H 2 = ^p- f --2H 2 f T , (19) 



2K 2 p + 6H 2 + f + 12H 2 f T 

24H 2 f TT - 2 - 2/ T ' 1 j 



where /t = df/dT, a prime denotes the derivative with respect to In a, and p, p are the total energy 
density and pressure, respectively. In an universe with only dust matter, p = p m = and p = p m . 
Obviously, if f(T) = const., f(T) theory reduces to the well-known ACDM model. 

In fact, f(T) theory was firstly used to drive inflation by Ferraro and Fiorini [24Ll25| . Later, Bengochea 
and Ferraro [2(|, as well as Linder [2l[, proposed to use f(T) theory to drive the current accelerated 
expansion of our universe without invo king dark energy. Very soon, f(T) theory attracted much attention 
in the community. We refer to e.g. |2rM}lj | for relevant works. 



III. CONSTRAINING f(T) THEORIES WITH THE VARYING 
GRAVITATIONAL "CONSTANT" 

A. \G e g/G e g\ for a general f(T) theory 



As is well known, a varying effective gravitational "constant" is one of the common features of most 
modified gravity theories [2j. Of course, as a modified gravity theory, f(T) theory is not an exception. 
In fact, the effective gravitational "constant" of f(T) theory has been derived in [30|, namely 



G c 



G 



N 



1 + ./7 



(21) 



Obviously, if f(T) is a linear function of T, namely f(T) = aT, the effective gravitational constant is 
just rescaled to be G e g = Gjv/(1 + ot), which is still constant in time. However, in general G e g is varying 
for any non-linear f(T). From Eq. (|2"Tj) . it is easy to get 



(22) 



Gog 




IttT 


G e fi 




1 + /t 



Using Eq. (|18p and the universal relation x = —(1 + z)H(dx/dz) for any function x, we have 

dE 2 



T = -(l + z)T H 



dz 



(23) 



where E 2 = T/T = H 2 /H 2 . Substituting f T , Jtt and Eq. ^ into Eq. @ and setting z = 0, the 
present value of |G e ff/G e ff| is on hand (note that the observational constrain given in Eq. ((T|) is obtained 
from the lunar laser ranging experiment, at redshift z = 0). It is easy to see that the Hubble constan Ho 
will appear in the final result of \G e s/G e g\ through T at redshift z = 0. According to [Hj], 



H Q = 100/ikm/s/Mpc = 1.02275 x l(T 10 /i yr" 1 



(24) 



Note that we use the units % = c = 1 throughout this work. It is easy to see that the Hubble constant Hq is 
a suitable measurement of the present value of \G e s/G e s\. Very recently, the SHOES (Supernovae and HO 
for the Equation of State) Hubble Space Telescope program [3J, [34[ released its latest model-independent 
measurement of Hubble constant, namely h = 0.738 ±0.024 (with only 3.3% uncertainty) (34[. Therefore, 
we adopt h — 0.738 in the followings. 

In the next subsections, we will consider two concrete f(T) theories, namely, f(T) = p(—T) n and 
f{T) = —P-T (l — e^ T °/ T ) , which are the most popular f(T) theories discussed extensively in the literature 
(see e.g. [IMEHM3)- 
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Description 


Best fit 


2cr right edge 


2a left edge 


2a top edge 


2a bottom edge 


(Q m o, n) 


(0.272, 0.04) 


(0.308, 0.04) 


(0.240, 0.04) 


(0.272, 0.26) 


(0.272, -0.29) 


\G e s/GeB | 


1.86849 


2.01124 


1.72109 


15.8991 


10.1177 



TABLE I: The present value of |G e ff/G e ff| (in units of 10 12 yr x ) with the best-fit parameters of [22J and the 
corresponding 2a edge for the case of f(T) = /i(— T) n . 

B. f(T) = ll (-T) n 

At first, we consider the case of f(T) = /.i(—T) n , where n and n are both constants. This is the simplest 
model, and has been considered in most papers on f(T) theory. Obviously, if n = 0, it reduces to ACDM 
model. Substituting f(T) — T) n into the modified Friedmann equation (|T9")) . one can easily find that 
/i is not an independent model parameter, namely (2lL l22l . l3lj | 



i - n 



M = 



2n 



^37- ( ~ To) 



where fi m o = K 2 p m o/(3i?o) ^ s the present fractional energy density of dust matter. So, we have 

m = ^%(-T ) 



2n- 1 



and then 



/t 



7t(l - flmo) £ 2(n-l) 



/tt 



1 - 2?1 

n(n-l)(l-r! m o)^ 2( "- 2) 



1 - 2n 



Tn 



fro 



Itto 



n(l - O m0 ) 
1 - 2n ' 

_ n(n - 1)(1 - n m0 ) 
(1 - 2n)T 



(25) 

(26) 

(27) 
(28) 



where E 2 = T/T a = H 2 /Hji Substituting f(T) = fi(-T) n and Eq. ([25]) into the modified Friedmann 
equation (fT9"]) , we find that [22l . l3l| 



£ 2 



n m0 (i + z) 3 + (i - o m0 )i; : 



2» 



(29) 



Obviously, if n = 0, it reduces to the one of ACDM model. Differentiating Eq. (|2"9"|) with respect to 
redshift z, we obtain 



dE 2 

dz 



3fi m0 (l + zf 



l-n(l-fi m o)£ 2 M) 



(30) 



G ff 




3n(n - 1)(1 - Cl m0 )Q, m0 H 


G e ff 




(1 - n)^ - n 2 n 2 m0 



Substituting Eqs. ([27]). (j2"gj) and ([301) into Ec l- <E3) and then Eq. (|2"2"]t. it is easy to see that at redshift 
z = 0, we have 



(31) 



Therefore, if the model parameters £l m o and n are given, we can correspondingly get the present value 
of \G cS /G cS \ from Eq. (JHH). Note that in [22|, this f(T) = fi{~T) n theory has been constrained by 
using the latest cosmological data, i.e., 557 Union2 type la supernovae (SNIa) dataset, baryon acoustic 
oscillation (BAO), and cosmic microwave background (CMB) data from WMAP7. The corresponding 2a 
results are given by [22| 



f^mo — 0.272^^ 



0.036 
032 > 



n = 0.04t ; 22 . 



(32) 



G 



0.04 



0.02 



S 



-0.02 



-0.04 




0.25 0.26 0.27 0.28 0.29 0.3 

&m0 



FIG. 1: The viable region in the Q m o — n parameter space which can satisfy the observational constraint on the 
varying gravitational "constant" given in Eq. ([TJ, as well as the latest cosmological data (SNIa+BAO+CMB), for 
the case of f(T) = ^{-T) n . See text for details. 



At first, we try to see whether the present value of \G e s/G e s\ with the best- fit parameters of [22| and the 
corresponding 2er edge can satisfy the observational constraint on the varying gravitational "constant" 
given in Eq. (JTJ, namely, \G/G\ < 1.3 x 10~ 12 yr -1 . From TableU it is easy to see that none of them can 
satisfy the observational constraint on the varying gravitational "constant" given in Eq. (pQ). Also from 
Table |H we find that for a fixed n, the smaller fi m o, the smaller \G e s/G e s\ is, whereas \G e g/G e s\ is also 
somewhat correlated with the extent of n deviating from zero. Next, we try to find the viable region in 
the fl m o ~ n parameter space which can satisfy the observational constraint on the varying gravitational 
"constant" given in Eq. (TTJ), as well as the latest cosmological data (SNIa+BAO+CMB). To this end, 
we scan the f2 m o — n parameter space within the 2er region of |22}. namely 0.240 < Cl m o < 0.308 and 
—0.29 < n < 0.26 (see Eq.(|3"2"])). and correspondingly calculate the present value of \G c s/G e s\ at every 
scanned (f2 m o, n ) point. The viable parameter region is determined by |G e ff/G e ff| < 1.3 x 10 _12 yr _1 , 
and we present it in Fig. [TJ which is approximately a trapezoid region. It is easy to see that the allowed 
f2 m o is still unchanged, i.e., 0.240 < f2 m o < 0.308. However, the allowed n is significantly shrunk from 
the original —0.29 < n < 0.26 of [22] to a very narrow range around zero, namely 

- 0.032 < n < 0.030. (33) 

This is the additional constraint on f(T) — /i(— T) n theory from the varying gravitational "constant". 
Clearly, it is very powerful in fact. 

C. /(T) = -fiT (l - e^ T °/ T ) 

In this subsection, we consider the case of f(T) = -fjtT (1 - e' 3To / T ), where fx and /3 are both constants. 
Obviously, when j3 — > we have f(T) — > fj,/3To = const., so it reduces to ACDM model. Substituting 
f(T) = — fj,T (l — e^ T °/ T ) into the modified Friedmann equation (TPS")) , one can easily find that \i is not 
an independent model parameter [2l|, [H, [3l[ , i.e., 

1 - fi m o , . 

M= l-(l-2/3) e ^ (34) 
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Description 


Best fit 


2cr right edge 


2a left edge 


2a top edge 


2a bottom edge 


(OnO, P) 


(0.272, -0.02) 


(0.308, -0.02) 


(0.238, -0.02) 


(0.272, 0.29) 


(0.272, -0.22) 


|G e ff/G e ff| 


0.919465 


0.988932 


0.842703 


9.74454 


13.8775 



TABLE II: The present value of \G eS /G eS \ (in units of 10" 



yr 



) with the best-fit parameters of [22| and the 



corresponding 2a edge for the case of /(T) = —pT I 1 — e 



J5TO/T 



It is easy to obtain 



/T = -M + M(l-^ 2 )e^ 2 



/to 



-/x + //(!- e* 5 



(35) 



TT 



VP . e 0/B 2 



Jtto = -=r e ■ 
Jo 



(36) 



where E 2 = T/T Q = H 2 /H 2 . On the other hand, substituting f(T) = -fiT (l 



Friedmann equation (jT9|) . we find that [22t l31j 

S 2 -r! m o(l + ^) 3 + M^ 2 



o/3/B 2 



) into the modified 



(37) 



If p — > 0, we have — > 1 — fi m o from Eq. (j34j) . and hence Eq. (|37|) reduces to the one of ACDM model. 
Differentiating Eq. (|3"T|) with respect to redshift z, we obtain 



dE 2 
dz 

Substituting Eqs. (JSHJ) , dHSJ) and 

z = 0, we have 



3n m0 (l + z) 2 



(38) 







G e ff 





1 + (j, [(1 - /3/£ 2 + 2/3 2 /^ 4 ) e^/ £2 - l] ' 

into Eq. (1231 and then Eq. (|22p . it is easy to see that at redshift 



(39) 



[1 - fi + fi(l - [1 - /x + fi(l - p + 2p 2 )eP] 



where fi have been given in Eq. Therefore, if the model parameters Sl m0 and j3 are given, we 

can correspondingly get the present value of \G e s/G e g\ from Eq. (j39j) . Note that in (22|, this /(T) = 
— /iT (l — eP T °l T ~^ theory has been constrained by using the latest cosmological data, i.e., 557 Union2 
type la supernovae (SNIa) dataset, baryon acoustic oscillation (BAO), and cosmic microwave background 
(CMB) data from WMAP7. The corresponding la results are given by 22] 



n 979+0.036 
u -^' z -0.034 i 



P 



-0 02+ ' 31 



(40) 



Again, wc firstly try to see whether the present value of \G e s/G e g\ with the best-fit parameters of [22j and 
the corresponding 2a edge can satisfy the observational constraint on the varying gravitational "constant" 
given in Eq. ([T]), namely, |G/G| < 1.3 x 10 _12 yr _1 . From Table [TTl it is easy to see that the first three 
points (whose p is close to zero) can satisfy the observational constraint on the varying gravitational 
"constant" given in Eq. ([1]), whereas the last two points (whose /3 is far away from zero) cannot. Also 
from Table HIl we find that for a fixed p, the smaller O m o, the smaller \G e g/G e s\ is, whereas \G e e/G e g\ is 
also somewhat correlated with the extent of P deviating from zero. Next, we try to find the viable region 
in the f2 m o — P parameter space which can satisfy the observational constraint on the varying gravitational 
"constant" given in Eq. (HJ, as well as the latest cosmological data (SNIa+BAO+CMB). To this end, 
we scan the f2 m o — P parameter space within the 2a region of [22j], namely 0.238 < Q m o < 0.308 and 
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0.04 



0.02 



03. 



-0.02 



-0.04 

0.24 0.25 0.26 0.27 0.28 0.29 0.3 

FIG. 2: The viable region in the Q m o — ft parameter space which can satisfy the observational constraint on the 
varying gravitational "constant" given in Eq. (TTJ, as well as the latest cosmological data (SNIa+BAO+CMB), for 
the case of f(T) = -^T (l - e f>TalT ^j . See text for details. 

—0.22 < p < 0.29 (see Eq. (|4T)|) ), and correspondingly calculate the present value of \G c s/G c e\ at every 
scanned (f2 m o, ft) point. The viable parameter region is determined by |Gcff/G c ff| < 1-3 x 10 _12 yr _1 , 
and we present it in Fig. ® which is approximately a trapezoid region. It is easy to see that the allowed 
ftmo is still unchanged, i.e., 0.238 < £l m o < 0.308. However, the allowed /3 is significantly shrunk from 
the original —0.22 < (3 < 0.29 of [22] to a very narrow range around zero, namely 

- 0.030 </?< 0.033 . (41) 

This is the additional constraint on f(T) = —pT (l — e^ T °/ T ) theory from the varying gravitational 
"constant". Clearly, it is very powerful in fact. 

IV. CONCLUSION 

As is well known, a varying effective gravitational "constant" is one of the common features of most 
modified gravity theories Q- Of course, as a modified gravity theory, /(T) theory is not an exception. 
Noting that the observational constraint on the varying gravitational "constant" is very tight, in the 
present work we tried to constrain f(T) theories with the varying gravitational "constant". We found 
that the allowed model parameter n or (3 has been significantly shrunk to a very narrow range around 
zero. The results are very impressive. 
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